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In Levin-Wen (LW) models, a wide class of exactly solvable discrete models, for two dimensional 
topological phases, it is relatively easy to describe single fluxon excitations, but not the charge 
and dyonic as well as many-fluxon excitations. To incorporate charged and dyonic excitations in 
(doubled) topological phases, a generalization of the LW models is proposed in this paper. We first 
enlarge the Hilbert space with adding a tail on one of the edges of each trivalent vertex, to describe its 
internal charge degrees of freedom. Then we study the full dyon spectrum of generalized LW models, 
including both quantum numbers and wave functions for dyonic quasiparticle excitations. The local 
operators associated with the dyonic excitations are shown to form the so-called Tube algebra, whose 
representations (modules) form the quantum double (categoric center) of the input data (unitary 
fusion category). In physically relevant cases, the input data is from a finite or quantum group 
(with braiding J?-matrices), we find that the elementary excitations (or dyon species), as well as any 
localized/isolated excited states, are characterized by three quantum numbers: charge, fluxon type, 
and twist. They provide a “complete basis” for many-body states in the enlarged Hilbert space. 
Concrete examples are presented and the relevance of our results to the electric-magnetic duality is 
discussed. 
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I. INTRODUCTION 

In recent years two-dimensional topological phases 
have received increasing attention from the science com¬ 
munity. These phases represent a novel class of quan¬ 
tum matter at zero temperature^, whose bulk proper- 
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ties are robust against weak interactions and disorders. 
They may be divided into two families: doubled (with 
time-reversal symmetry, or TRS, preserved), and chiral 
(with TRS broken). Chiral phases were first discovered 
in integer and fractional quantum Hall (IQH and FQH) 
liquids. Mathematically, their effective low-energy de¬ 
scription is given by Chern-Simons gauge theory or, more 
generally, topological quantum field theory (TQFT)^!^. 
Doubled topological phases include topological insulators 
and some states in quantum spin liquids. Either chiral 
or doubled phases may be exploited to do fault-tolerant 
(or topological) quantum computing^^—. 

The (chiral) Chern-Simons theories are formulated in 
the continuum and have no lattice counterpart. On 
the other hand, doubled topological phases do admit a 
discrete description. The first such formulation in the 
physics literature was the Kitaev’s toric code model^. 
(In the mathematical literature, a discrete version of 
TQFT had been constructed a bit earlier by Turaev and 
Viro^, which by now is known to describe certain doubled 
phases.) About ten years ago. Levin and Wen (LW)^ con¬ 
structed a wide class of discrete models on a trivalent lat¬ 
tice/graph, with an exactly solvable Hamiltonian, for two 
dimensional doubled topological phases. The model is 
now believed to be a discretized version of doubled Chern- 
Simons theoryiS, which is mathematically the same as 
the Turaev-Viro TQFT— The original motivation of 
the LW model was to generate ground states that exhibit 
the phenomenon of string-net condensationi^ as a physi¬ 
cal mechanism for topological phases. The ground states 
in this model can be viewed as the fixed-point states of 
some renormalization group flowi^, which look the same 
at all length scales and thus have no local degrees of free¬ 
dom. Like Kitaev’s toric code model^, we expect that the 
subspace of degenerate ground states in the LW model 
can be used as a fault-tolerant code for quantum compu¬ 
tation. 

Two of us have studied, in a previous joint paper with 
another author—, the ground state degeneracy (GSD) of 
the LW model on a (discretized) closed oriented surface 
M. Usually in TQFT the GSD is examined as a topolog¬ 
ical invariant of the 3-manifold x in the LW 

Hamiltonian approach, our computation of the GSD be¬ 
came accessible to physicists. In this paper, we attack the 
problem of solving the full spectrum of quasiparticle exci¬ 
tations in the LW models with the input data being a uni¬ 
tary fusion category (see below for details). This problem 
is of significance for further interdisciplinary study of the 
models in physics, mathematics and quantum computa¬ 
tion codes. It is generally believed that the quasiparti¬ 
cle excitation species are related to the quantum dou¬ 
ble that classifies the degenerate ground states. Several 
proposals about excitation spectrum in the LW models 
have been made in the literatureSii2ri£. In this paper, 
we will present a new approach to understanding the full 
elementary excitation spectrum of the LW models, that 
addresses both the quantum numbers and corresponding 
states (or wave functions) explicitly for all quasiparticle 


species. In particular we want to accommodate the needs 
in physics and quantum computation codes for concrete 
and explicit expressions to play with. 

Several new developments feature our analysis. Usu¬ 
ally for a single (pure) fluxon, it is easy to character¬ 
ize/specify their quantum numbers. See, for examplej^S. 
However, the fusion of two or more (pure) fluxons, gener¬ 
ally leads to the appearance of charge quantum numbers. 
(Some examples are shown later in Sec. VH.) Namely the 
set of single fluxon species are not closed under fusion! 
So how to represent all dyon (charge-fluxon composite) 
species in the LW model presented a challenge. Our way 
to solve this problem is to enlarge the original Hilbert 
space of the LW model. We explicitly introduce internal 
charge degrees of freedom (d.o.f.) at each trivalent vertex 
by adding a tail on one of its edges. This has greatly facil¬ 
itated the treatment of fusion outcomes. In this way, the 
LW model is actually generalized, with the underlying 
graph(s) involving univalent vertices. The second impor¬ 
tant development is that we have identified the operator 
algebra for the local operators, that can be used to gener¬ 
ate all quasiparticle excitations, to be the Tube algebra 
constructed by Ocneanu22r^. Using the Tube algebra, 
the relationship of quasiparticle species to the irreducible 
representations (simple modules) of the quantum double 
becomes relatively easy to establish, and the above men¬ 
tioned complicated situations for fusion of non-abelian 
anyons become easier to handle. It has been showii2^ that 
a half-braiding in the quantum double (or center) cate¬ 
gory corresponds to an irreducible representation (simple 
module) of the Tube algebra. This not only enables us 
to define the string operators, but also to account for 
charged as well as dyonic excitations. Our present anal¬ 
ysis has clarified and emphasized the importance of sup¬ 
plementing the twist, as quantum number in addition to 
the usual charge and fluxon-type, to the characterization 
of quasiparticle species in generalized LW models. In¬ 
deed, generally there may exist quasiparticle excitations 
which have the same charge and fluxon-type but have 
different twists and, therefore, should be counted as dif¬ 
ferent species. 

A similar operator algebra approach for quasiparti¬ 
cles in topological phases has been proposed by Lan and 
Wei>i^. The Q-algebra they independently invented is 
a presentation of the Tube algebra. In their approach, 
they added an extra (charge) index at each vertex for its 
internal charge d.o.f., while we want to add a (charge) 
tail at one of the links attached to a vertex. This dif¬ 
ference makes the two formulations have complimentary 
technical advantages and disadvantages. 

Because of the interdisciplinary interests in the LW 
models, we have tried to adapt our presentation in this 
paper to audience with different backgrounds. Of course, 
the basic audience in our mind is physicists, and we have 
tried hard to make the presentation accessible to physi¬ 
cists. However, whenever a reference of the terminology 
or of the idea can be made to the mathematical litera¬ 
ture, we will do it to help readers of mathematical back- 
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ground. (Readers with physics background can safely 
skip these mathematical remarks without harming their 
further reading.) The last section is also devoted to the 
relationship between our approach and TQFTs. 

We will use some terminology in category theory lan¬ 
guage for convenience of physicists, because this language 
could be used widely in future physics, just like group 
theory has become the language of contemporary physics. 
Condensed matter physicists do not need to be worried. 
Whoever has learned angular momentum or crystal group 
theory in quantum mechanics is familiar with at least 
one fusion category, which is nothing but the category 
formed by all finite-dimensional unitary representations 
of the rotation group or its discrete subgroups in three- 
dimensional Euclidean space! The decomposition of the 
(tensor) product of two irreducible representations into 
a (direct) sum of irreducible representations just gives to 
the fusion algebra, with the non-negative integral coef¬ 
ficients in the direct-sum decomposation as the fusion 
(rule) coefficients. The 6j-symbols are well-known in 
group theory. So the fusion category is a straightforward 
generalization of the representation theory of groups (or 
group algebras) to more complicated algebras (more pre¬ 
cisely, weak Hopf algebras). Up to now only fusion cate¬ 
gories associated with a finite group or a quantum group 
appear in the literature of condensed matter physics. 

The paper is organized as follows. In Section II we 
briefly review the LW models and set up our notations. 
In Section III we review the topological symmetry of the 
ground states, i.e. the invariance under Pachner muta¬ 
tions of the (spatial) graph. Then we begin our study of 
excited states by introducing an extension of the Hilbert 
space of the LW models, as well as the extended Hamil¬ 
tonian in Sec. IV. We devote Sec. V to the central issue 
of the paper, i.e. the study of elementary (quasiparti¬ 
cle) excitations, using local operators preserving topo¬ 
logical symmetry, which is shown to form the tube al¬ 
gebra. Minimal projection operators and simple mod¬ 
ules (irreducible representations) of the tube algebra are 
introduced. A dyon species is identified with an irre¬ 
ducible representation (simple module) of the Tube alge¬ 
bra, and fusion of all dyon species gives rise to the quan¬ 
tum double (or the categoric center) of the input fusion 
category. String operators are generalized to dyon-pair 
creation, hopping operators etc, and their properties that 
are related to important observables, such as twist and 
5'-matrix, are studied using graphic calculus. Next two 
sections. Sec. VI and Sec/ VH, are devoted to study¬ 
ing excitation spectrum and emergent braiding statistics 
from the above-established set-up. In Sec. VHI, with 
possible a pplications in anyon condensation, we exam¬ 
ine the particular case - the braided LW models - with 
the input fusion category equipped with an i?-matrix. 
Physically these models are actually a generalized gauge 
theory with gauge “group” being a finite or quantum 
group. In Sec. IX, we present several examples, includ¬ 
ing cases with input data from an abelian group, from a 
non-abelian group S 3 , from the Kitaev’s quantum dou- 



FIG. 1: A configuration of string types on a directed triva- 
lent graph. The configuration (b) is treated the same as (a), 
with some of the directions of some edges reversed and the 
corresponding labels j conjugated j*. 


ble model as well as from a modular category, such as 
the double semion model and double Fibonnaci model. 
Sec. X addresses the electric-magnetic duality between 
two particular LW models with two seeming different sets 
of input data involving the same finite group. In Sec. XI 
we elaborate the relation to topological quantum field 
theory for the convenience of readers with mathematical 
background. The final section (Sec. XH) is devoted to 
conclusions and discussions, in particular, we emphasize 
that the local operators we have defined in this paper 
and their algebra may be useful even if the Hamiltonian 
is deformed away from the LW Hamiltonian. 


II. LEVIN-WEN MODELS 

Let us briefly review the Levin-Wen models. The input 
data to define the model, i.e., to specify the Hilbert space 
and the Hamiltonian, is the unitary fusion category C. 
More specifically, we will use the tensor description of C 
in terms of 6j-symbols. 

The model is defined on a trivalent graph embedded to 
a closed oriented surface. The Hilbert space is spanned 
by the degrees of freedom on edges. See Fig. [TJ For 
each edge, we assign a label j (called string type), which 
runs over a finite set of integers L = {j = 0,1,...,V}. 
The Hilbert space is spanned by all configurations of the 
labels on edges. Each label j has a “conjugate” j*, which 
is also an integer and satisfies j** = j. If we reverse the 
direction of one edge and replace the label j by j* on 
this edge, we require the state to be the same. See Fig. 
[H There is unique “trivial” label j = 0 satisfying 0* = 0. 

To specify the Hamiltonian, we introduce the structure 
on string types as follows. 

A fusion rule on L is a function NiLxLxL^N 
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such that for a, b,c,d€ L, 

NL = N^aO = Sab, (1) 

N°b = Sab^, (2) 

= (3) 

x^L x^L 

A fusion rule is multiplicity-free if £ {0,1} for all 
a,6, c £ L. We restrict to the multiplicity-free case 
throughout this paper unless specified. We define Sabc ■= 
N^b which has the symmetric properties: Sabc = Sbca and 
Sabc = Sc»b»a*- We say a triple {a,b,c) is admissible if 
Sabc — f- 

Given a fusion rule on L, a quantum dimension is a 
map d : L —>■ R such that da* = da and 

^ ^ dc^^abc* — dadfo. (4) 

c&L 


In particular, do = 1. Let aj = sgn(dj) which takes 
values of ±1 for each label j, and require: 


When q = i, the equality implies that must varnish 

unless SmiqSk*ip = 1- By using the tetrahedral symmetry, 
one arrives at Eq. ([7]). Here Vj = -^/dj is a choice of a 
square root of the quantum dimension. The number Vj 
is either real or pure imaginary, depending on the sign 
aj = sgn(dj). 

Depending on how the square root is taken, Vj is de¬ 
termined up to a sign. We fix the sign as follows. From 
the conditions in Eq. ([S]), we have (G'^Q^jVjVk)^ = Sijk- It 

is possible to fix the sign of Vj such that Gg^jVjVfe = Sijk- 
We define 


'-^ooi 


( 9 ) 


In particular, vq = I because do = 1 (from Eq. (|4])) 
and thus Gqqq = 1 from Eq. ([6]). Indeed, we can verify 
v^ = dj directly from the orthogonality condition in Eq. 
([6|) together with do = 1. The definition in Eq. m also 
implies 


G 


ijk 

Okj 


^j^k — ^ijk ; 


( 10 ) 


CXiG-jG-k — I, if Sijk — 1- 


( 5 ) 


Given a fusion rule and a quantum dimension on L, we 
may define 6j-symbols, often denoted as G. A tetrahe¬ 
dral symmetric unitary 6j-symbol is a map G : L® —>■ C 
satisfying the following conditions: 


'-'kin ~ '-'nk-l* ~ '-'ijn* ~ GmGn Orj.j,. 


mlq 

Lrr, 


E l ^mlq 

n kp* n'^pk* n 


mns 
I* m*i 


^js n _ ^jip ^riq 
^Ikr* — ^q*kr*^mls*'> 

_ Siq r r 


( 6 ) 


The data {dj, <5^^, can be derived from the rep¬ 

resentation theory of a group, or more generally a quan¬ 
tum group. (More generally, such a set of data is from a 
unitary fusion category.) For instance, we take the labels 
j to be the irreducible representations of a finite group. 
The trivial label 0 is the trivial representation. The fu¬ 
sion rule tells whether the tensor product ji 0 j 2 0 js 
contains the trivial representation or not. The number 
Gj is the Frobenius-Schur indicator telling if the represen¬ 
tation j is real or complex, or pseudoreal, dj = ajdim(j) 
the dimension dim(j) of the corresponding representa¬ 
tion space multiplied by the Frobenius-Schur indicator 
Gj. The number G^")™ the (symmetrized) Racah 6j sym¬ 
bol for the group. In this example, the LW model can be 
mapped to the Kitaev’s quantum double model. 

One important property of the 6j-symbols is that 

(7) 

To prove this, one can rewrite the orthogonality condition 

by 


= S,qSmuSk*v, 


which can be proved by the pentagon identity 
boGo),jG^c(°Ggo\° = and the orthogonality 

, r-njk _ J_r. 

'^j'-'Qkj'-'k* 0 j ~ 

There are two types of local operators, Qy defined at 
vertices v and (indexed by the label s = 0,1,..., N) 
at plaquettes p. Let us first define the operator Qy. On 
a trivalent graph, Qy acts on the labels of three edges 
incoming to the vertex v. We define the action of Qy on 
the basis vector with ji, j2, ja by 



( 11 ) 


where the tensor Sj-^j^j^ equals either 1 or 0, which deter¬ 
mines whether the triple {ji, j2, js} is “allowed” to meet 
at the vertex. Since iJjijats i® symmetric under permuta¬ 
tions of the three labels: Sj-^j-^h = Sj^j^j^ = Sj-^j^j^, the 
ordering in this triple {ji,j2,j3} is not important. 

The operator B® acts on the boundary edges of the 
plaquette p, and has the matrix elements on a triangle 
plaquette. 



=ViiVj2Vj3VjjVj^Vj^G 


isilia 


^3^32 3i ^3 g33 32 

•^iii^* ^i^ia* 


( 12 ) 


The same rule applies when the plaquette p is a quadran¬ 
gle, a pentagon, or a hexagon and so on. Note that the 
matrix is nondiagonal only on the labels of the boundary 
edges (i.e., ji, j 2 , and js on the above graph). 

The operators B® have the properties 


b;^ = b;\ 

B;B; = Y,SrsfBi, 

t 


( 8 ) 


(13) 

(14) 
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Both can be verified by using conditions ([6]). 

The Hamiltonian of the model is (here D = d|) 

H = — Qy — Bp, Bp = — dgHp (15) 

V p S 

where the sum run over vertices v and plaquettes p of the 
trivalent graph. 

The main property of Qy and Bp is that they 
are mutually-commuting projection operators: (1) 

[Qv,Qv'] = 0 = [Bp, Bp'], [Qy,Bp\ = 0; (2) and QyQy = 
Qy and BpBp = Bp. Thus the Hamiltonian is exactly 
soluble. The elementary energy eigenstates are given 
by common eigenvectors of all these projections. The 
ground states satishes Qy = Bp = 1 for all v,p, while the 
excited states violate these constraints for some plaque¬ 
ttes or vertices. 

In particular, if {d, S, G} arises from the representation 
theory of groups or quantum groups, we have Syst* = 
Ssrf- Then the Bp's commute with each other, 

[b;,,b;^] = o (16) 

which can be verified by the conditions in ([6]) when pi 
and p 2 are the two nearest neighboring plaquettes, and 
by Eq. (ITdll together with Jrst* = Ssrt^ when pi = P 2 - 

III. TOPOLOGICAL SYMMETRY FOR 
GROUND STATES 

To characterize the topological phases, we study the 
topological observables. Examples include the topologi¬ 
cal degeneracy of ground states. Behind them the topo¬ 
logical symmetry plays an important role: topological 
observables are those invariant under mutations of the 
spatial graph. In continuum theory, they are observables 
invariant under the smooth deformation of the space-time 
manifold. In the following we analyze the mutation sym¬ 
metry for the ground states. 

Let us begin with any two arbitrary trivalent graphs 
and discretizing the same surface (e.g., a torus). 
It is known that they can be mutated to each other by a 
composition of the Pachner moves^i: 






FIG. 2: A mutation two graphs that discretize the same man¬ 
ifold. The left one is mutated to the middle one by a compo¬ 
sition of fi moves, and the middle one is mutated to the right 
one by a /s move. 


To the elementary moves /i, / 2 , and fs, we associate 
linear maps between the corresponding Hilbert spaces as 
follows: 


Ti : 


T. : 


n : 





J3 


J '3 ’ 


J4:J5J6 






Vd 


^4 363^ 


32 


( 20 ) 

)' 

( 21 ) 

( 22 ) 


Note that since we can reverse any edge by conjugat¬ 
ing the corresponding label, the above formulas do not 
depend on the edge directions. 

Between the Hilbert spaces and on any 

two graphs, there is a mutation transformation by a 
composition of these elementary maps. In particular. 
Bp = D~^ '^sBp is a special example. In fact, on the 
particular triangle plaquette p as in (1221) . we can verify 
Bp=s/ = E2T3. 

The ground states have properties: 


1. The mutations are unitary in the ground-state sub¬ 
space. 


2. The ground states are invariant under mutations. 
By unitarity Oi = we mean {cj)\ Oi \cj)') = 

{(p'\ O 2 Eor example, = T 3 can be verified by 



T 2 


h 


Vji Ujs Uj-g 


_ Is le ^^312 i * _ 
yCp 3i3s3t 


■h 


T. 


3 ^ 


-03 ' 


Ayis \ (23) 


to 2 


using condition ([6]), XliX 


See Eig.[2]for instance. 

We can associate two different Hilbert spaces to P*-^^ 
and respectively, as described in the previous sec¬ 
tion. Denote by T-LB) the Hilbert space on P^^^ and 
on r(2). 


IV. EXTENSION OF THE MODEL 

To study the spectrum of Levin-Wen models, we will 
extend the Hilbert space. An elementary excitation 1-0) 
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(a) 


(b) 


FIG. 3: Extension of the Hilbert space by extra tails. 



lattice in Fig. |3(a)[ There are three edges connected to 
each vertex. To each vertex, we associate an open edge 
called a tail and attach it to on e of the three connected 
edges. For example, in Fig. |3(b)] each vertex carries a tail. 
To define the Hilbert space it does not matter which one 
of the three edges neighboring to the vertex we choose to 
attach the tail. 

The Hilbert space is spanned by the string types j £ I 
on all the edges of the tailed graph. With the tails, there 
are two more d.o.f. around each vertex. For example, in 
Fig. |4(b)l these extra d.o.f. are labeled by k and q near 


is actually presented by two vertices in Fig. |3(b)[ Around 
each new vertex in Fig. |3(b)[ we require the fusion rule 
Sijk = 1 for the three neighboring edges labeled by i,j, 
and k connecting to the vertex. For example, at the left 
upper corner of the plaquette in Fig. |4(b)[ we require 

^jlhkl = 1 S-nd Skij*q* = 1. 

The Hamiltonian has two terms: 

H = -J2Qv-J2^p- (24) 

V p 

The first term is 


each vertex. Each vertex of the spatial graph in Fig. 3(a) 




FIG. 4: Around each vertex, two extra d.o.f. are needed: q is 
assigned to the tail, and k to the line connecting the vertex 
and the tail. 



(25) 


supports two types of quasiparticles, charge at vertex v 
if Qy \tp) = 0, and fluxon at plaquette p, if Bp {ip) = 0. 
We extend the local Hilbert space at v to support distin¬ 
guished charges and internal degrees of freedom between 
them. 

Let us start with a trivalent graph, e.g., the hexagonal 


The second term is 

= ( 26 ) 

s 

with 


J 





E 


Vj3 Vj, Yj, Yj, Vj- Vj/ Vj. Yf^ Yf^ V^-. X 


Vfc4Vfe5Vfc' ^ 



( 27 ) 
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The operator is a straightforward generalization of 
the Levin-Wen operator in Eq. (na). The two tails 
labeled by and <75 are viewed as two external legs of the 
plaquette. With qi = 0 = ( 72 , the plaquette is effectively 
treated as having 8 boundary vertices. Acting on this 
effective plaquette by 73® defined in Eq. (EH) , we arrive 
at 16 of v’s and 8 of 67 -symbols in the Eq. (EZl). 

If we restrict to the Q = 1 subspace, we recover the 
traditional Levin-Wen model. Take Fig. |4(b)| for an ex¬ 
ample, in the subspace with <7 = 0 fixed for all tails, we 
have ki = j 2 = ^ 2 ,^ 4 = Js = h, h = ke = Iq, etc. 
Hence, Q = 1 subspace in Fig. 4(b) is identified with the 
usual Q = 1 subspace in Fig. 4(a) In such a case, H® 
becomes the usual 73®. 

The model is exactly solvable because the local terms 
in the Hamiltonian (12411 are mutually commuting projec¬ 
tion operators. 


V. ELEMENTARY EXCITATIONS 

We study elementary excitations by algebra of local 
operators preserving topological symmetry. Excitations 
support quasipartciles, which are identified with the irre¬ 
ducible representations of the algebra and are classified 
by the quantum double category. 

Elementary excitations are mutual eigenvectors of all 
Q’s and H’s. In particular, the ground states are Q = 
1 and B = 1 eigenvectors. Since the Q = 1 subspace 
recovers traditional Q = 1 Levin-Wen Hilbert subspace, 
the ground states in the extended model are exactly the 
same as in the traditional Levin-Wen model. 

Elementary excitations support local quasiaprticles. If 
an excitation is a = 0 eigenvector, we say there is 
a charge quasiparticle living at u, which is identified by a 
nontrivial tail label. On the other hand, if Sp = 0 (with 
Qv = i around the plaquette) we say there is a fluxon 
quasiparticle at p. We call a generic quasiparticle a dyon 
- a composite of charge and fluxon. 


out tails. The ground states have mutation symmetry as 
discussed in Sec IHIl In the following we explore the topo¬ 
logical symmetry in excitations. The corresponding topo¬ 
logical observables under such symmetry form the “tube 
algebra”. This enable us to classify the elementary exci¬ 
tations. In particular, the good quantum numbers can be 
identified by the irreducible representations of the tube 
algebra, which are formulated by the quantum double of 
the input category. 

Excitations have less symmetry. In the presence of a 
nontrivial quasiparticle at the triangle plaquette in Eg. 
dUD and ( 1 ^ . excited states are not invariant under Tj 
and Ts, because T 2 T 3 = Bp=^ = 0 . 

In the generic case, the tail label — called a charge — 
within a plaquette is nontrivial. States are not invariant 
under T 2 and T 3 . But they still have Ti symmetry. An¬ 
other symmetry with respect to the tail is to move a tail 
along the plaquette boundary. Define 



(28) 


which has a similar form to Ti. The tail <71 can freely 
move along the plaquette boundary, as long as it does 
not cross another tail, e.g., the one labeled by <72 above. 


A. Topological symmetry and tube algebra of 
observables 

As discussed above, when all tail labels are trivial, the 
extended Hilbert space becomes the traditional one with- 


We are interested in topological observables, i.e., local 
operators that preserve Ti and T 4 transformations. Con¬ 
sider elementary excitations with at most one quasipar¬ 
ticle at the plaquette for simplicity. Hence we consider 
only one tail for simplicity as follows. Define 


















Vil Vj3 Vj, Vj, Vfc Vj J Vj/ Vj/ Vj/ Vj/ Vj/ Vfc' X 



= E 
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^sj[k'* ^sr^3[* ^ 



( 29 ) 


This operator has a graphical presentation of fusing a 
string labeled by s along the plaquette boundary, per¬ 
formed by T transformations (up to some normalization 
factor), see Fig. O 





FIG. 5: Graphical interpretation of Bgjg/g and Bp. (a) Bgsq'u 
attaches a string and fuses it along the plaquette boundary 
by Ti and T 3 . (b). With g = g' = 0, B^jqi^ is reduced to 
BoaOs = Bp. (c) With w = 0, 0 performs a rotation of a tail 
along the plaquette boundary. 

The topological observables preserving Ti and T 4 are 
linear combinations of Bpkqt- Denote such a generic op¬ 
erator by 


Which fuses a loop to the plaquette boundary. See Fig. 

EKb). 

Quasiparticles in elementary excitations are identified 
by projection operators 11 • If = 11 where 11 is minimal. 
Here minimal means that if H is a combination of projec¬ 
tion operators H = Hi -|- n 2 then Hi or II 2 is zero. Each 
projection H projects onto states with a specific quasipar¬ 
ticle at p called a dyon. A ground state is a rip^p — 1 
eigenstate; we say such a p has a trivial dyon which we 
identify with the projection Bp. Dyons identified by all 
other projections II Bp carry higher energy, because 
Bp is a special minimal projection and Bp ■ II = 0. 

Consider dyons with charge q at the tail fixed. Oper¬ 
ators on such states form a subalgebra Aq with elements 
^ = J2kt^<ikqtBqkqt- If we express a projection in Aq by 

n, = 'Y2 AqktBqkqt, (33) 

kt 

then Aq-Iiq = H, implies 

n,,, = d,d, ^ (34) 


— ^ ^ ^pkqt^pkqt-) 
pkqt 


(30) 


where the summation runs over p, k, q, t with 6 pqt* = 1 = 
^kqf (otherwise Bpkqt = 0). These operators satisfy the 
multiplication rule x ■ y = z where z is given by: 


Zpkqt = dkdt ^ xipqrypmisGZ4 ( 31 ) 


mnlrs 


which can be verified in graphical presentation in Fig. 
[SKa). The operators in Eq. (l30l) equipped such a multi¬ 
plication rule form the tube algebra A. 

The local operator Bp in the Hamiltonian is a special 
element in A'. 


Bp — dgHosOs- 


(32) 


Each minimal projection H^ identifies a dyon at p and 
is in one-to-one correspondence with an irreducible rep¬ 
resentation of Aq. 


B. Quantum double theory of dyons 

The dyons with fixed charge q are identified with irre¬ 
ducible representations of Aq. However, the dyon species 
are identified with the irreducible representations of tube 
algebra A. A dyon with fixed charge q does not form 
a species itself, because a generic topological observable 
Bqkq't transforms the dyon with charge q to other dyon(s) 
with charge q'. A dyon species is identified with a set 
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of dyons that are invariant under A, i.e., with an irre¬ 
ducible representation of A. Such irreducible representa¬ 
tions form quantum double category, with each represen¬ 
tation identified with a quantum double label. Hence the 
natural algebraic theory of dyons is the quantum double 
category theory. 

In the rest of this subsection we expand on the ideas of 
the previous paragraph. The irreducible representations 
of tube algebra form a quantum double category. The 
key structure in the latter is the half braiding. A half¬ 
braiding tensor z satisfies the naturality condition 


E j j sl* f^s*pm^rn*tr* 

'^r^s^lnqrZpmls^nr*t ^jn*t ^q*n*k 


(35) 

A minimal solution to this equation is associate with 
a quantum double label J. Each J is one-to-one 
corresponding to an irreducible representation of tube 
algebrai^^. Denote each minimal solution by . 

Quantum double labels classify dyon species. Each 
dyon species may carry different charges q just like each 
spin may carry different magnetic components. In this 
case, we say these dyons belong to the same dyon species, 
denoted by a quantum double label J. A dyon species 
is identified by a minimal central projection in A. Here 
“central” means it commutes with all topological observ¬ 
ables in A. It is a sum 


n^ = En,^ (36) 

<? 


where q runs over all dyons that belong to J. 

Each dyon that belongs to species J has the projection 
Hg arising from by 


H*^ 

H*^ 


dfcdt J 

, ^qkqf 
Ug 


(37) 


There may be several projections Hp, Hg, ..., arise from 
the same J, with p ^ q. 

In general, each J may carry multiple copies of a charge 
q, denoted by index a in H*^ = J2q a ^q,a- Throughout 
this paper for simplicity, we assume each q appears at 
most once in all J. 


labeled by the trivial string type q = Q. We draw the 
dotted line to present the trivial label 0 for convenience. 
Fix an edge e, and consider a state jT) with no charge at 
either of the two vertices of e. For such an edge we can 
define a creation operator. For example, in the follow¬ 
ing diagram, we create a pair of dyons across the middle 
vertical edge labeled by ■ Define creation operator by 



The resulting state I'k) has charge p at the 

lower left vertex, and charge q* at the upper left vertex. 
The two plaquettes are occupied by a pair of dyons. 

The generated state is normalized to: 


(^1 p^^J;pg*t^J;pg* ^ dpdg/dj(^' jd-) . (39) 


Given a ground state, enables us to explicitly 

write down the elementary excitation wavefunction. 

Now we develop the dyon string operator, which cre¬ 
ates a pair of dyons at the end of the string. First, we 
create two dyon-pairs across two edges respectively, e.g., 
the two labeled by jg and by jg as follows. 


C. Dyon string operator 


In this subsection we define dyon-pair creation and an¬ 
nihilation operators. In the Q = 1 subspace, all tails are 
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The resulting state has the first pair of dyons occupying 
the left and the middle plaquette, and the second pair 
occupying the middle and the right plaquette. Next we 
annihilate the two charges in the middle plaquette. By Ti 
followed by a sequence of T 4 moves, we can move the tail 
q' to be at the same edge with g'*. Then we annihilate 
the charge by 



The desired string operator is the composition of the 
following: we first create two dyon-pairs as in Eq. ioi), 
second we move the two tails in the middle plaquette to 
be at the same edge, and third we annihilate the two 
charges as in Eq. m- This procedure defines a string 
operator across two edges. We can repeat this procedure 
to define a string operator along a longer string. 

The process to annihilate charge pairs q and g* and 
sum over g at the middle plaquette is called the con¬ 
traction of charges. After the contraction, no nontrivial 
quasiparticle is left at the middle plaquette. The string 
operator defined is path independent: two string opera¬ 
tors along two isotopic paths result in the same final state 
if the final position of the dyon is the same and there is 
no non-trivial quasiparticle in the area enclosed by the 
two paths. 

We use ribbon strings to represent creation and string 
operators. In Fig. [ 6 l^a), the ribbon string represents the 
creation by of a dyon pair across an edge e. There 

are charges p and q* at the two ends of this string. In Fig. 
[SJb) the dotted line presents the contraction of charges at 
the middle plaquette. Here the charge contraction con¬ 
nects two strings to a new one. The two string operators 


in Fig. El^b) and Fig. [ 6 l^c) are equal, illustrating the path 
independence of the string operator. 


D. Twist and modular S matrix 

In this subsection, we analyze some topological observ¬ 
ables in terms of the string operators to characterize the 
dyon species. 

We need a special choice of to write down the cre¬ 
ation operator (or equivalently the elementary excitation 
wavefunction). This amounts to picking a specific repre¬ 
sentation in which the state transforms under the tube 
algebra A. On the other hand, the topological observ¬ 
ables that characterize topological properties of dyon ex¬ 
citations do not depend on specific choices of z'’. In 
the following, we explore the topological properties using 
n*^, which are uniquely determined by the 6 j-symbols 
and does not depend choices of 

The simplest invariant is obtained by contraction of 
the charges at the two ends of a string, leading to a loop 
operator. If no nontrivial quasiparticle exists in the area 
enclosed by the loop, the closed string operator gives a 
multiple of the identity matrix: 



where dj is called the quantum dimension of J, defined 

by 

= (43) 

qeJ 

The next topological observable is the twist: define the 
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by the string operators 



Since the string operators are path independent, it is not 
important where on the spatial graph we put the string. 
The crossing matters, which indicates the path order of 
creation and charge contraction operators. Therefore it is 
safe to draw only the ribbon strings without mentioning 
the underlying spatial graph. 




FIG. 7: A dyon twist is presented by a twist of the ribbon 
string. 


FIG. 6: (a) A ribbon presents a pair of dyons created at 
two plaquettes across edge e. (b) Two dyons are created. 
The dotted line presents the contraction of the charges at the 
middle plaquette. The resulting state is dyon pair state with 
two dyons at the two ends of the ribbon string, (c). Creation 
along two isotopic strings result in the same dyon pair state if 
no nontrivial quasiparticle exists in the area enclosed by the 
two strings. The two string operators in (b) and (c) are the 
same. 


Another important topological observable is the mod¬ 
ular S matrix, as defined as follows. First, we create 
a dyon-pair. Second, we add a closed string operator 
around on end point of the previously created dyon-pair. 
Third, we contract the ends of the dyon-pair. See Fig. [S] 
If no nontrivial quasiparticle exists in the area enclosed 
both strings, the final operator is a multiple of the iden¬ 
tity matrix. Presented in terms of the ribbon strings, 
they are 


twist by 




SjK^- 


( 47 ) 


e = ( 44 ) 

9 


It commutes with all dyon projections, and hence is a 
good quantum number of a dyon state. For states with 
dyon identified by 11;^, the eigenvalue is solved to be 


0j = 


1 


EnL 


(45) 


The S matrix characterizes the holonomy effect of 
winding dyon J around K, or equivalently, exchanging 
J and K twice. 

The S matrix turns out to be independent of choices 
of z'^. It evaluates to be 



The twist 9j is the same for all dyons in the same species 
J, even with different charges q. 

This scalar is a 17(1) number which we identify with dj. 
The definition in Eq. (|4^ has a graphical presentation 
of a self rotation of the tail. See Fig. EKc). Hence 6 
is identified with the dyon’s statistical spin s via 6 = 
exp(27ris). 

If we apply the twist in Fig. HKb), we obtain a string 
operator as in Fig [T] Twisting either dyon in the mid¬ 
dle plaquette of Fig. HKb) before the charge contraction 
leads to the string operators in FiglTJa) and FiglT^b) sep¬ 
arately. Therefore the twist Oj of a dyon can be detected 


Sometimes it is useful to define T matrix by 

Tjk = SjK0j- ( 49 ) 

Modular matrices S and T are believed to characterize 
the quantum double category, and contain all kinematic 
information of the dyon species. 

E. Fusion and hopping operators 

Dyon species are closed under fusion. Here the fu¬ 
sion process is described as follows: when two pairs of 
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FIG. 8: Three steps to evaluate S matrix: (a) create a pair 
of A'-dyons; (b) wound around a K-Ajon by J-dyon closed 
string operator; (c) contract the two ends of the A'-string. 


dyons of species J and K are created on the same two 
plaquettes, the resulting state is a linear combination of 
ones obtained from the creation of one pair of dyons. If 
the dyon-pair state L appears in this linear combination, 
then we define 5jkl* = 1 and 5jkl* = 0 otherwise. 


Next we consider another way to describe the fusion 
process that results in elementary three-dyon states (on 
the sphere). In terms of string operators, we create three 
pairs of dyons of species J,K and L, with three dyons 
at one end of each pair meeting at the same plaquette. 
If we annihilate these three dyons, the resulting state is 
the zero vector or a nonzero three-dyon state. See Fig. 
Ha). We define the fusion rule by 5jkl = 1 if we obtain 
a nonzero three-dyon state and Sjkl = 0 otherwise. (In 
general, there may be more than one fusion channels, 
but throughout the paper we consider the multiplicity 
free cases for simplicity.) 


Let us consider the later fusion process in more detail. 


To do this we define the fusion of charge by 



where the charge q 2 is moved upward and fused with qi, 
resulting a linear combination of charge q[ state. 

The fusion process is illustrated in Fig. Ha) and de¬ 
scribed as follows. First, we create three pairs of dyons 
by IT'^’^-^ , and with summation over p, g 

and r. Second, we annihilate the three charges at the 
middle plaquette by composition of Ti moves and the 
above charge fusion given in Eq. (I50p . Then we apply 
the projection Qy. Finally, we apply Bp at the plaquette 
to annihilate the fluxon. If Sjkl = 1 , the nonzero result¬ 
ing state is graphically presented by a riboon three-valent 
tree structure. See Fig. Hb). 

The fusion rule is completely determined by the S ma¬ 
trix (known as Verlinde formula3^): 


Sjkl 


1 ^ SjnSknSln 
^ 5'oAr 


(51) 


We end this subsection by describing the hopping oper¬ 
ator: First, we create two dyon pairs, both of species J, 
with two particular charges k and k* are at the same pla¬ 
quette. See Fig. [TUI Second we annihilate charge using 
Eq. (IdTI) and fluxon using Bp. 

The hopping operator is equal to the composition of 
with the charge and fluxon annihilation. We 
present this equality graphically by 





(52) 


VI. EXCITATION SPECTRUM 

In Sec. IV El we showed how to generate a three-dyon 
excitation from a ground state. In this section we study 
the full spectrum for all excitations. 
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FIG. 9: (a). Create three pair of dyons, and fuse three dyons 
at the middle plaquette by annihilating charges and fluxons. 
The dotted line presents charge annihilation and dashed line 
presents fluxon annihilation, (b). If Sjkl = 1, the fusion 
process results in a elementary three-dyon excitation. 



FIG. 10: Hopping Operator: the dash circle represents fluxon 
annihilation. 

We consider models on the sphere only. There is only 
one ground statei^. We claim all excitations can be gen¬ 
erated from the ground by string operators and annihi¬ 
lation operators. 

To reveal the structure of the spectrum, it is convenient 
to consider unit cells of one vertex (including one charge) 
and one plaquette. For simplicity, let us hrst consider 
a simplified situation: each plaquette contains at least 


one tail. At each plaquette, we only consider nontrivial 
charge at a particular tail. This amounts to enforcing 
all other tails at each plaquette to be labeled by trivial 
charge. Then each plaquette together with the tail form 
a unit cell that supports exactly one dyon. See Fig. fTlT al 
for example. 




FIG. 11: (a). Each plaquette has one tail that may take 
nontrivial charge, (b) Unitary Ti moves mutate the graph to 
a tree-like graph, with P the total number of plaquettes. 

By applying Ti moves, one can always mutate the 
graph to have a tree-like graph as in Fig. [TlT bl. All pla¬ 
quettes are transformed to P — 1 bubbles except the last 
one. (The outside region forms one plaquette on sphere.) 
These bubbles are the new unit cells that support ex¬ 
actly one dyon. All (good quantum numbers of) dyons 
at bubbles are preserved during the mutation. 

A typical excitation can be generated by P — 1 pairs 
of dyons across the bubbles. Let us denote the dyons 
inside bubbles by {Jp,qp}p<p-i. But the P — 1 dyons 
outside the bubbles form a huge multiplicity. Similar 
analysis based on fusion process in Sec. IV El implies that 
this multiplicity can be organized by tree structure: 

Jiqi J292 J393 Jp-2qp-2 Jp-iqp-i 



The fusion channels that occur outside region are diag¬ 
onalized by {Ke}i<e<p- 3 - These K degrees of freedom 
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describe how the P — 1 dyons are fused into linear com¬ 
bination of Jpqp dyons at the outside plaquette. 

Therefore the basis for excitations is 

{\{Jpqp}l<p<P,{Ke}l<e<P-3) I 

(54) 

^JlJ2K*^Kp-3Jp-iJp ^KeJe+2^e + l 

e 

In the models arising from modular tensor categories, 
each quantum double label is a pair ij. See Sec. IIX Cl 
The basis is simplified as 


{ I 5 Qp}3<p<P^ {^e^e}l<e<P—3) I 

P-4 

p e 



VII. EMERGENT BRAIDING STATISTIGS 


The basis (l54l) allows us to calculate the fractional ex¬ 
change statistics of dyons. The transformation of de¬ 
generate 7V-dyon states under the exchange of any two 
dyons can be computed using the hopping operators we 
have developed in Sec. IV El They form a representation 
of the Braid group B]\j, because of the path independence 
of the hopping operators. 

Consider 7V-dyon excitation states, with N dyons la¬ 
beled by {Jp( 7 p}i<p<Ar at N fixed unit cells (plaquette 
together with a tail). The braiding matrix is computed 
in the iV-dyon excitations have basis 

{|{-^"e}l<e< 7 V- 3 ) 

(56) 

\SJiJ2K*Skn-3Jn-iJn 

e 

Although they dyons may be different of different species, 
the braiding matrices form a presentation of the braid 
group Bn. The braiding matrices are independent of 
charges qp. 

In models arising from modular tensor category, any 
three dyon states has the basis 

I*iji,<?i;*2j2,g2;i3j3,<?3) • (57) 

The braiding matrices cti and a 2 are diagonal matrices 
given by diagonal of and 

Consider the doubled Fibonacci model for example. 
Consider the four-fluxon states on a sphere. Each fluxon 
is either a pure fluxon labeled by (J = rr, g = 0), or 
carrying a charge r, labeled by (tt). The four-fluxon 
states have a basis 


TT TT TT TT 



where K = 1 , r, r, tt. The four dots at the top label four 
fluxons, which may or may not carry a charge r. These 
charges do not affect the braiding matrices, and are thus 
not presented in the basis. For simplicity, we choose all 
four fluxons to be (J = tt, g = 0), and the computation 
is within the usual Levin-Wen Hilbert space. 

If we exchange two fluxons in the counterclockwise di¬ 
rection by the hopping operators, we obtain the braiding 
matrices in the above basis 

/ I 0 0 0 \ 

0 0 0 

ai - (T 3 — , 

0 0 e 5 0 

V 0 0 0 I / 

/ 02 e-^03/2 ^ ^ 

e“'^0^/2 e“^02 0 e^0^/2 

^ 0^/2 0 gv 02 e ~^ 0^/2 ’ 

\ 0 0^/2 e~^ 0^/2 02 J 

(59) 

where 0 = . ai exchanges the fluxon I and 2, 0-2 

exchanges 2 and 3, and 0-3 exchanges 3 and 4. They 
generate the representation of the braid group B 4 . 

The four eigenvalues of ui and <73 are verihed to be 
B^- <72 can be obtained by a basis transformation in 

terms of 6j symbols. 


VIII. BRAIDED MODELS 


Many example models are equipped with B matrix, 
include the models arising from representations of finite 
groups and quantum groups. The presence of an B- 
matrix simplihes the disruption of the operators in the 
model. In this section we analyze in detail how to charac¬ 
terize the dyons by three good quantum numbers: charge, 
fluxon, and twist. We study in more detail the creation, 
annihilation, and string operators in this situation. 

Let {d,(5, G} be the data discussed in Sec. [ITl The B 
matrix is a map B : ^ C that satisfies hexagon equa¬ 

tions: 


E l picad* pe pyabg _ pd p^acd* p/ 

^g^be* g^gc^ce" f ^ac^be" f ^bc 

9 

E l /-le^bdrye y-ie* a g pd y-ie*bd jjf 

^9^cag -^ad^bcf -^ac^acf ^ab' 


(60) 

(61) 


The data {d, <5, G, B} is a tensor description of a uni¬ 
tary braided category. Examples include the models aris¬ 
ing from representations of finite groups and quantum 
groups. 
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A. Good quantum numbers of dyons 


basis transformations among them: 


1. Charge 


Recall the definition of a charge: an excited state |'0) 
has a charge at vertex n, if IV') = 0, namely, if the 
tail label associated to the vertex v is not the trivial 
0. We say |'0) carries a charge q^. More precisely, define 



It commutes with the Hamiltonian (j24p . and thus q is 
good quantum number of {ip). In particular, = Qv 

projects onto trivial charge. 

Another good quantum number in the charge excita¬ 
tions is related to the topological spin of the charge. To 
construct it, we examine how a tail charge is associated 
to a vertex. There are different choices to associate a tail 
to vertex to specify the Hilbert space. In this section, 
there is no canonical choice better than the others. For 
example, if the three edges incoming into one vertex are 
labeled by i, j, and Z, there six possible ways to associate 
a tail labeled by q: 



All of six choices are equivalently good. These six ways 
specify a basis of six different Hilbert spaces. Define the 



The transformation v moves the tail between the left 
side and the right side on the same edge, while /i moves 
a tail to another edge. Both moves are in the clockwise 
direction. 

In the absence of fluxon at plaquette, the twist defined 
in Eq. SMh can be reinterpreted as 

0 ^ = h'fih'fih'fj,, ( 66 ) 

which take a tail in the counterclockwise direction around 
the vertex and finally back to the same position. This 
process realizes “self-rotation” of the charge. 

The eigenvalue for charge q is 


0 , = ( 67 ) 

which is a C/(l) number that depends only on q. 


2. Fluxon 

Suppose there is a nontrivial charge at plaqueete p. 
An excited state {ip) has a pure fluxon at plaquette p, 
if Bplpj) = 0. To identify fluxons in the presence of 
nontrivial charges on the tail inside p, we extend by 
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2^ Vjl Vi2Vi3Vi4Vi5Vj6VfciVfc3Vfc4Vfe,VjJ Vj-Vj-Vj/Vj/Vj/ X 


y^^334,32 ^^^ 4^414 ^^53e^5 /^^&3\3& ^ 
Vfc'jVfe^Vfe^Vfe^G^,^-,^,. s*k'^j’^* S*j!^j'/ s*j'^k’^ s'fcsie* s*-'-'* 


^6^1 


/ 0^2 D -?2 ^ 92 ^ 2^2 0^2 \ ^ 94 ^ 5^4 

s*A: 5^_72* <?i^iy y ^ 2^2 ^* 32 ^ 2 * ^232 J s*k'^j'^ s*j'^k'^ 



fej 

k' 

-*• qi ^5 

32 

35' 

k^ 




( 68 ) 


The operator S® is a straightforward generalization of the 
Levin-Wen operator in Eq. (112^ . It can be obtained as 
follows. We first apply the basis transformations ly on all 
tails that point into the plaquette p, i.e., the tails labeled 
by <71 and q 2 in the above example. The resulting graph 
contains the plaquette with four tails pointing outwards: 



Now the Levin-Wen operator is well-defined on this 
new plaquette, which is treated as having 10 vertices on 
the boundary with 10 external lines labeled by h, I 2 , 
..., /e, as well as < 71 , q 2 , ( 74 , and q^. After applying 
Bp, we move back the two tails qi and <72 towards right 
by the inverse transformations v~^. Loosely speaking. 
Bp = BpVlV 2 ^ where vi ^2 is the half twist on the 

tail < 71 , 2 . This results in Eq. (l 6 ^ . 

The formula in Eq. (1551) can be read as follows. The 
effective plaquette in (l69l) is treated as having 10 bound¬ 
ary vertices. According to the definition of B® in Eq. 
( 1121 ), we arrive at 20 of u’s and 10 of 6 j-symbols in the 
Eq. (l27l) . The four copies of the half twist v gives rise to 
four of R tensors (or their complex conjugates B), in the 
two brackets, for two tails <71 and q 2 respectively. 

All local operators B® and Ql are mutually commut¬ 
ing with any other. According to the interpretation in 
Eq. ((551) . this is a direct consequence of the property 


that Qy and Bp are mutually commuting projection op¬ 
erators. The operator B® can be recovered from B® by 

around p 2^)- 

In the following we use B® to identify fluxons at p. The 
ocal operators B® forms the fusion algebra 

&p&p=Y,5rst>Bl. (70) 

t 

with multiplication obeying the fusion rule. 

In the following we derive a set of orthonormal pro¬ 
jection operators in the fusion algebra to identify parti¬ 
cle species of the fluxons. The braided model equipped 
with R matrix has Sijk = Sjik- The algebra (1701) is now 
abelian, and it is uniquely determines a NxN matrix 
Ay, called the fusion characters, satisfying 

Ay = Xf (71) 

Xfxf = J2s^3k-X^X^ (72) 

k 

^Ay^ = 5^.B,^Ay^ = y,. (73) 

3 A 

The matrix Ay is unique up to the relabeling of A = 
0,1,..., A — I. Ay can be viewed as normalized one¬ 
dimensional irreducible representations of the fusion al¬ 
gebra, as observed in Eq (l72)l . The factor A^ on the 
RHS of Eq (172]) normalizes Ay to satisfy Eq. (1731) . 

The matrix Ay determines a set of orthonormal pro¬ 
jections operators at p: 

:= A/AyBy (74) 





















17 


satisfying 

n^n^=6A,Bn^, ^n^ = l. (75) 

A 

These projection operators identify the particle species 
A of the fluxons at p. Each projects onto the states 
with A-type fluxon at p. There is a special fluxon type, 
denoted by A = 0, corresponding to quantum dimensions 
by Xj = dj/VO. For A = 0, rip = Bp projects onto 
states without any nontirvial fluxon at p, and thus we 
say A = 0 is the trivial type. Each A comes with a 
conjugate A* such that = X^, and we say A*-type 
fluxon is the antiparticle of A-type fluxon. 


We define twist as follows. Take the unit cell of a vertex 
V and the plaquette p, and consider a dyon that carries 
charge q and fluxon A. The idea of the twist is to wind 
the charge q around the fluxon A at plaquette. 


3. Twist 


In additional to the charge q and the fluxon A, there 
is another good quantum number, the twist, which arises 
from exchange between q and A. For example, in the 
Z 2 gauge theory (toric code model), there are four types 
of elementary quasiparticles: the trivial one 1, the Z 2 
charge e, the Z 2 flux m, and the charge-flux composite 
em. The twist of em is —1 because the wave function 
acquires the Arharonov-Bohm phase —1 by exchanging 
e and m twice(or equivalently, by winding e around m 
once), which renders em the fermionic statistics. 


Start with the plaquette as in Fig. |4(b)[ and consider 
the dyon living at the unit cell of the vertex of qi and 
the plaquette. Suppose it carries the charge qi and the 
fluxon A. We shall construct the twist operator to move 
the charge qi around the plaquette. This operation is 
similar t o Bf, w hich creates a Wilson loop labeled by s, 
see Fig. 12(a) However, the twist differ by moving qi 
along an open path labeled by gi, with the open end 
treated as the new tail, see Fig. |12(b)[ 



FIG. 12: (a) Bp creates a Wilson loop labeled by s; (b) the 
twist of the charge qi moves qi around the plaquette. 


Define the twist by 
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'?1^2.?3* 'Ji^sie* 'JlJe-?!* ^*^1^2* \ '?2^1 (7*fcifc2* *?2-?2y *?i^‘4y5* ^iJs 



• (76) 


When the fluxon is trivial at plaquette p, 0vp becomes 
the special twist 0„ of the pure charge at v in Eq. (1551) . 
Therefore 0„p takes the dyon around the entire unit cell 
vp. 

It turns out that twist Qyp is commutes with both 
and Bp and yields a good quantum number. A dyon 
in an excitation \tjjj) is characterized by the three good 
quantum numbers, the charge q, the fluxon A, and the 
twist 0j, if 

Ql \^j) = IV'j), 

np |iAj) = IV-j), (77) 

0vp IV'j) = li’j) ■ 


Two dyons that carry the same charge q and the same 
fluxon A may have different twists, because the twist 
measures more information than the Wilson loop, see 
Fig. 12(a) and Fig. |12(b)| The Wilson loop acting of 
the dyons is completely determined by the test charge s 
and the fluxon A of the dyon. It creates a pair of charge s 
and s *, winds s around the fluxon A and then annihilate 
s and s*. The entire process gives rise to the Arhoranov- 
Bohm phase X^/Xq according to Eq. ([74]). The twist, 
on the other hand, exchanges the charge q and the fluxon 
A twice, and measures more information about the states 
that is characterized by the twist 6 j. 

Although our discussion is restricted at a particular 
plaquette in Fig. |4(b)[ the definition of the twist 0„p is 
valid on any unit cell of a vertex v and a plaquette p 
(with V on the boundary of p). 


B. Dyon-pair state 



FIG. 13: (a) A circle with two two-valent vertices on the 
sphere, with two vertices vi,V 2 and two plaquettes pi,P 2 - (b) 
Circle with two tails. 


good quantum numbers qualitatively as well as quanti¬ 
tatively. 

Start with a circle on sphere with two two-valent ver¬ 
tices, see Fig. 13(a) The Hilbert space is spanned by 
the for degrees of freedom on the circle with two tails, 
denoted by p*, q, x, and y, see Fig. |13(b)[ We use p* for 
the future convenience. 

We divide the space into two unit cells vipi and V 2 P 2 - 
The excitations are classified by the two dyons living at 
these two unit cells. The two dyons are always paired. 

We shall explore the following properties of these dyon- 
pair states in the section: 


1. No single (nontrivial) dyon exists on a sphere. 

2. If there is no (nontrivial) fluxon, the charges at vi 
and V 2 are opposite to each other. 

3. Two fluxons at pi and p 2 are opposite to each other. 

4. Two dyons have the same twist. 


In this subsection, we study the lowest excitation on 
sphere, the dyon-pair states. We shall study the three 


The charge projections at the two vertices vi,V 2 are 

(78) 


Ql\=Sp 


Qll = s. 
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The fluxon projections at two plaquettes pi , p 2 are 

( 79 ) 

S 


with being 



and 



(80) 


(81) 


Let us fix the fluxon A at pi, and consider the dyon- 
pair states with charges p* at vi and q at V 2 fixed. Such 
dyon-pair states are Up^ Ql^ = 1 eigenstates. The 
total number of distinguished dyon pair states is then 

=trin^,Ql,QV 

= ^ d,dyX^GU/,Glfpy ,, ( 82 ) 

xys 


where in the second equality we used property that 

UX — A . . 

-ttpyllpy — Up*Xy*- 

In non-abelian models, i.e., with |djj > 1 for some 
string type j, the charges p* and q may not match exactly 
as p = q. However, there is some selection rule to pair p* 
and q, depending on the fluxon A. 

We prove the properties as follows. When the fluxon 
at Pi is trivial; H = 0, the dyon excitations are rip^^ = 1 
eigenstates. If we fix the charge p at Vi and g at U 2 , from 
Eq. ([82|), the total number of possible states is 


tr(n^^ =^^d:,dydsG‘i4AGPJ 


D ^ 

xys 


(83) 


=5, 




Next we prove property 3. By definition of R matrix, 
the expression in the bracket in Eq. (l80l) can be expressed 
as 



We also rewrite the expression in the bracket in Eq. m 
similarly, followed by the substitution = R%g, and 
Rj^g, = Ry»g. After substituting the formula in the 
brackets, we find 

( 85 ) 

which implies 

^ti=^P 2 - (86) 

If the dyon at vipi has fluxon A, then the fluxon of the 
dyon at V 2 P 2 must be the its anti-fluxon A*. 

It can be also proved that 


0 


vipi 


= 0 


V 2 P 2 ■ 



(87) 


( 88 ) 


The two dyons carry the same twists. 

Elementary excitations are dyon-pair states. In each 
pair the two dyons have the same twists, opposite fluxons. 
The charges of two dyons may not match exactly, but 
satisfy some constraint that depends on the fluxon. 

An elementary excitation is given by an simultaneous 
eigenvector tjj oi , Q 2 J ^ 1^2 s-nd 0 i _2 with eigenvalues 
being the quantum numbers p*,q,A^0. They are given 
by the half-braiding tensors 



'^x^yZpt^ljtySpp'Sqqi , 


(89) 


for some quantum double label J, where p* and q are 
fixed charges that two dyons carry. 

Each quantum double label J is parameterized by 
fluxon type A and the twist 6 . 

In a special case, the ground state is 


X 




(90) 


C. Creation, annihilation, and string operators 


where in the second equality we used Eqs. © and O- 
If the charge p* at vi is fixed, there exists one (and only 
one) dyon excitations with charge charge at V 2 being p. 
This proves the property 2. In the particular case when 
p = 0, then the only allowed state is the ground state, 
and hence property 1 is proved. 


The properties of dyon-pair states analyzed in the pre¬ 
vious subsection hold on a generic graph. To see this, 
here we shall consider the creation and annihilation op¬ 
erators of dyons and string operators, in the setting of an 
i?-matrix. This will enable us to generate all elementary 
excitations from a ground state. 
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1. Charge string operators 

We first study the pure charge case. Recall the ground 
states are the simultaneous = 1 and Bp = 1 eigen¬ 
states. In the Q = 1 subspace, all tails are labeled by 
the trivial string type q = 0. As above, we draw the 


dotted line to present the trivial label 0 for convenience. 
In the following we give an explicit formula for creation 
operator that in terms of the R matrix. 

In the Q = 1 subspace the creation operator that cre¬ 
ates a pair of charges at the two ending vertices of an 
edge e by 



Here e denote the left boundary edge of the plaquette. 
The resulting state is an pair of charges, q at the top 
vertex, and q* at the bottom vertex. 

The operator is normalized as follows: 

(V^l IV') = d, (^1^). (92) 


If the two tails on the edge e are not on the same side, 
W§ is defined up to a basis transformation /i or iz acting 
on Eq. dSH). For example. 


J 



In cases where i?-matrix exists, each charge qo itself 
forms a dyon species with trivial fluxon type. The corre¬ 
sponding half braiding tensor is 

^pkqT ~ ^p,qo^q,qo^qok- ( 94 ) 


Eq. (IMll is a special case of Eq. (l38l) . In general, without 
a R matrix, a charge does not form a dyon species. 

If there are nontrivial charge already present inside the 
plaquette, is defined by 





32 32 ^2 ^2 9 ^2 9|<?2 

3W1Q2 



(95) 
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which creates two charges: q at the upper left vertex and 
q* at the lower left vertex. The operator in Eq. (l95l) 
is unitary. One sees that when qi = 0 and q 2 = 0, Eq. 
m is recovered. 

The operator can be used to recover some of the 
previously defined operators. In particular, the charge fu¬ 
sion operator define in Eq. (1501) is equal to Q 2 
Also, the special case when qi = q^ = q in Eq. (1051) then 
Q 1 Q 2 J2q' dq'Wg' is equal to charge annihilation opera¬ 
tor define in Eq. 62). 

The above shows that the hopping operator defined in 
Sec. IV El can be express in terms of the R matrix. Also, 
the operators in this subsection allow us to define the 
string operators completely in terms of R matrix. 


2. Fluxon string operator 


In this subsection, we study the string operators for 
pure fluxons that carry no charge. We restrict to the 
Q = 1 subspace. 

We define the creation operator on an edge e by 






(96) 


It is diagonal in the matrix form. Only two plaquettes 
are shown, assuming the rest of the graph is unaffected. 
The definition holds for arbitrary shaped plaquettes. 

The operator generate a fluxon-pair state from 
any ground state |<I>), with fluxon A* on pi and A on p2, 
where pi is plaquette left to the edge e and p 2 right to e: 

1$) = 5 a , |d>) 

n^,W,^\<^>) = 5B,oWf\^). (97) 

These properties can be proved using the conditions (l6|) 
on 6j-symbols. 

The definition of does not depend on the direction 
of the edge e. In fact, if we reverse the direction of e, 
je in Eq. (|M)l is replaced by j*. implies 

where e and e~^ are the same edge with 
opposite direction. Both and W^i create the same 
fluxon pairs across the edge, see Eig. [Ml 

Erom Eq. (IMl), is the identity operator when 
A = 0, as creating a trivial fluxon pair does nothing. 
The hermitian of creates a conjugate pair of fluxons 
because XJ^* = X^: 


Wf = 


(98) 


In general (non-abelian case, i.e., with |dj| > I for 
some j in the input data), even a pure fluxon carries 
charges. The operator is a special case of a generic 
fluxon creation operator tlM="4;00 quantum double 
label J = A and with trivial charges at both ends. 




FIG. 14: Fluxon-pair state |4>) generated from a ground 
state 1$). The creation operator does not depend on the edge 
direction. The fluxon-pair state |<l?) in (a) is the same as 
Wf, 1$) in (b). 



(a) (b) (c) 


A 



(d) 


FIG. 15: (a) Three neighboring plaquettes around a trivalent 
vertex, (b) Create two fluxon pairs across the edge 62 and 
€ 3 . (c) Annihilate fluxons at pi by (d) The final fluxon- 
pair state in (c) is equal to that obtained by directly creating 
a fluxon pair across edge ei. This implies is path 

independent, and thus is a hopping operator of fluxon A at 
Pi- 


In the following we show how to annihilate and hop 
fluxons in the absence of charge at the plaquette. 

Let us start with a ground state |d>), and consider a 
trivalent vertex and its three neighboring plaquettes po, 
Pi and P 2 , see Fig. &). In the following, we suppress 
by W/^ for i = l,2, 3. 

In Fig. ITKI bL creates a A*-A fluxon pair at po 
and pi, while creates a A*-A fluxon pair at pi and 
P2- Now Pi is occupied by two fluxons, A from W^, and 
A* from The resulting state may no longer be an 

eigenstate of certain n® , because A and A* may couple to 
more than one types of fluxons. The operator |d)) 

can be decomposed by the orthonormal projections . 
The operator projects onto the state |<i)) 

with only i?-fluxon at pi. 

Particularly, kills any nontrivial fluxon at pi. In 
the above example, projects onto a fluxon-pair state, 
with A* at Po and A at p 2 - In this killing process, 71°^^ 
plays the role of annihilation operator. The annihilation 
can occur only if the two fluxons at pi are antiparticles 
of each other. 

The above process is also a hopping process, in which 
the hopping operator moves the fluxon A from pi 

to p 2 - In this process, a A-fluxon is created at p 2 while 
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(a) 

FIG. 16: 


a ^-fluxon is annihilated at pi. 

The hopping operator must satisfy some topological 
property: hopping along two homotopic paths (without 
any nontrivial quasiparticle enclosed by the two paths) 
leads to the same final state. Consider again the above 
example. We apply the hopping operator to the 

fluxon pair state |$), and obtain a fluxon pair state. 
The path independence requires 

(99) 

around any trivalent vertex. This property can be veri¬ 
fied by using the conditions ([ 6 ]) on 6 j-symbols. 

The hopping operators induce a string operator that 
creates a pair of fluxons far apart. We choose a path 
along plaquettes pi,p 2 , ■ ■ ■, and Pn+i, going across edges 
ei, 62 ,..., and e„, as illustrated below: 


Vi 

V2 



ei 

£2 


This is a string consisting of plaquettes. First we create a 
fluxon pair on the neighboring plaquettes across ei, with 
A*-fluxon at Pi and ^-fluxon at p 2 . Then we move the 
A-fluxon to Pn by a sequence of hopping operators, and 
the final state is 

( 101 ) 

The two fluxons are at the starting plaquette pi and the 
ending plaquette Pn+i of the string. The string oper¬ 
ator in Eq. (I101|l only depends on the two ends of the 
string because of the path independence of the hopping 
operator. 

IX. EXAMPLES 

There are many examples of input data for the models 
considered in this paper. Including examples related to 
the representations of finite groups, the group algebra of 


finite groups, and the representations of g-deformed uni¬ 
versal enveloping algebra of Lie algebras. In this section, 
we discuss some typical examples of these three classes. 

A. From finite group representations 

Given a finite group H, Levin-Wen models admit two 
different types of input data: from representations of iJ, 
with labels identified as irreducible representations; and 
from the group itself, with labels identified as group ele¬ 
ments. We call the former the Repn model and the latter 
the VecH model. 

In this subsection we consider several examples of mod¬ 
els arising from representations of a finite group H. To 
this end we now discuss a few general features in this 
context. 

The models are based on a tensor description of the 
representation category Repn of H. String types j are 
(representatives of) irreducible representations {pj,Vj). 
Quantum dimensions dj = ajdim(Vj) are equal to the 
dimension of the representation space, multiplied by the 
Frobenius-Schur indicator aj, which is 1 if the represen¬ 
tation j is real or complex, and — 1 if pseudo-real. 

The fluxons are classified by the conjugacy classes. 
Since the number of conjugacy classes is equal to the 
number of irreducible representations, the number of 
fluxons is equal to the number of charges, as expected 
from the analysis in previous section. 

Let {C^}a be the set of conjugacy classes of H indexed 
by labels A. The fusion characters are just the usual 
characters Xj(^) for H (up to normalization factors): 

where \H\ is the order of iJ, and \C^\ is the cardinality 
of C^. Note that XqXj = a^dim)!^-) = d^. The or¬ 
thogonality relations CSl) for X^ result from those for 
character functions. 

The quantum double labels are classified by pairs 
(A, p), where A labels a conjugacy class of H, and p is an 
irreducible representation of the centralizer Za = {ff S 
H\ghA = hAg}- Here Ha is a arbitrary representative 
element in but fixed once and for all. 


1. Abelian group 

Consider an abelian group H. All irreducible repre¬ 
sentations are 1-dimensional, and hence dj = 1. The 
6 j-symbol is given by 

(103) 

Each group element is itself a conjugacy class, so the 
quantum double labels are pairs {g, p) of group elements 
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and irreducible representations of H. Each dyon is a 
charge-fluxon composite. 

For example, let i? = Z^v, the quantum double charges 
are (g, fi) for g, /r = 0,1,..., iV — 1 and the z tensors are 

exp{2Trig/N)6pjt-Sjgt *, (104) 

where Spjt* = liip + j — t = 0 mod N and 0 otherwise. 


The nontrivial R matrix is i ?22 = ~1- 

There are three conjugacy classes, labeled by A = 
0,1,2, with \C^\ = 1,2,3 respectively. The fluxons are 
classified by the three conjugacy classes, with the char¬ 
acter table presented in Table HI 


TABLE I: Character table ol H — S 3 . 


2 . Reps^ model. 

Consider the model arising from representations of 83 . 
The string types are the three irreducible representations 
of the symmetry group S 3 , denoted by L = {0,1, 2}. All 
labels are self-dual, i.e., j* = j. The fusion rules are 
given by i5ooo = <5oii = 1^022 = <^122 = 1^222 = 1 - 

The quantum dimension dj is the dimension of the 
representation space Vj : do = <^1 = 1 and ^2 = 2. The 
independent nonzero symmetrized 6 j-symbols are 

1 

71 ’ 

1 

(105) 

All other nonzero 6 j-symbols are obtained through the 
tetrahedral symmetry in Eq. ®. 


^000 ~ -‘-5^111 — -^7^222 — “7^5^011 ~ -*-5^222 — 


1 

71 ’ 


^022 _ 1 ^022 _ 1 ^022 _ 1 ^122 „ 1 ^122 _ 

'-'022 — 2’^122 — 2’^222 — 2 ’ ^122 — 2 ’ ^222 — 


Xi(A) 

^A=0 ^A=l ^A=2 

Xi=o 

111 

Xi=i 

11-1 

Xi=2 

2-10 


There are 8 quantum double labels. Indeed, the cen¬ 
tralizers for the three conjugacy classes are Za=o = 83 , 
Za=i — Z 3 , and Za =2 — ^ 2 - In total there are 8 ir¬ 
reducible representations of Za. We denote 8 quantum 
double labels by J = I, 2,..., 8 . 

We present dyon-pairs graphically by a string 
pt=:4=»q^ with fluxon A. All distinguished dyon-pair 
states and the corresponding twists are enumerated in 
Table HD 

The properties developed in previous section can be 
verified, e.g., the total number of dyon-pair states for 
fixed A,p,q obey the counting formula in Eq. (15^ . 


0 

II 

t-H 

II 

0 

i 

0 

02 = 1 

2^=4^2 03 = 1 

A = 1 

^ A^l ^ ^ A^l 

A^l ^ A=1 ^ ^ 

2«=^^2 05 = exp(2|i) 

2«=^^2 06 = exp(-2|l) 

(M 

II 

^^ A^2 A^2 

A=2 ^ A=2 ^ 

^ ^ A^2 ^^ A^2 

A—2 A—2 ^ 



TABLE II: 17 dyon-pair states in Reps^ model. 

The explicit wavefunction for each dyon-pair is specified by the half-braiding tensors z: 


^pjqt — ^pP^q.p 


Zpjqt — ^p,25q,2 




^10 0^ 


^ 0 1 0 \ 

Z'pjqt — Spp5q,Q 

0 1 0 

+ SppSqA 

10 0 


[ 00 - 1 ) 


lo 0 i/ 


— ^i6pfi6qA6j^3St^3 + -^f(5pp(5q^o<5j,3<^t,3; 


jt 
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^pjqt ~ '^P,2'^i3:2 


^pjqt ~ ^p,2^q,2 


^pjqt ~ ^pfiSqfi 


^pjqt ~ '^P.1'^9,1 


/ 0 0 1 

0 0-1 
_ 2i7r ITT 2i7r 

ye 3 e 3 e 3 

/ 0 0 1 

0 0-1 

2i7r _ ITT _ 2i7r 

e 3 e 3 e 3 


ji 


jt 


10 0 
0 —1 0 I + 

ji 


\0 0 0 

^010 

— 1 0 0 I + <^p,2<^q,2 

V 0 0 0/,. 


^0 0 1 

0 0 1 I + SpfiSq^2Sj^3St^3 + Sp^2SqfiSj^3St^3j 

Vi 1 0/,. 

^ 0 0 1 \ 

0 01 + iSp,l^q,2Sj,3^t,3 + i^p,2Sq,lSj,3^t,3- 

V-1 -1 0/.. 


B. Kitaev’s quantum double model 

Here we consider the Levin-Wen model arising from a 
finite group H itself. Set the string types to be group 
elements: I = with h* = h~^. Set d/i = 1, for 

all h G H and Sghk = 1 if ghk = 1 and 0 otherwise. Set 

^kln = S^jmSklm-Sjkn^Sinl- (106) 

Hence = 1. 

The model is identified with Kitaev’s quantum double 
model on the dual triangulation graph. The local opera¬ 
tors form a quantum double algebra D{H) of H. 

Let 6^ = i3fc-igfe,fe-i,fc~ig for q,k G H. The tube alge¬ 
bra has the multiplication rule: 

bPb‘^, = (107) 

which recovers D{H). 

The dyons in elementary excitations are determined 
by solutions of Eq. (IMl) . Fix q at the tail, set = 
nq,fc-yfc-iq for k G Zq = {t G H\tq = qt}. The 
equation tt^ = Emez, ^li^'L-ik the solution tt^ = 

Xa{k) given by the character of all irreducible rep¬ 
resentations a oi Zq. 

The dyon species are identified by a pair (A, a) where 
A is a label of a conjugacy class of H and a G 
Irrep(Zq) for a representative element q in The mod¬ 
ular matrices are 

5'(A,a),(B,/3) = Xa(h)xjl(g), (108) 

^ ^ gGA,heB 

gh=hg 

T(A,a),{B,p) = for any g G A. (109) 

Both RepH model and the quantum double model have 
excitations classified by the same quantum double labels 


(A, a). This reveals an electric-magnetic duality: the 
former support quasiparticles of charges at vertices and 
fluxons at plaquettes while the latter support charges at 
plaquettes and fluxons at vertices. 


C. From modular category 


For a braided model with input data {d, S, G, i?}, as in 
Sec. IVIlIi defines the S matrix 

Sab = J2d,R:,Rt,. ( 110 ) 

c 

If S matrix is invertible, the input data is a tensor de¬ 
scription of a unitary modular category. The quantum 
double classification is quite simple in this case. 

The quantum double labels are pairs denoted by ij, 
with quantum dimension 


= (Ill) 


The fluxon type of ij is j. In particular, pure fluxons are 

jj- 

The half-braiding tensors are 


dadbRtkR^Cg'^Glii^, . ( 112 ) 


ab 


The S matrix for the quantum double labels are 

= S.kWi. ( 113 ) 

and the twist is 

e-^=Rl,/Rl,. (114) 

Modular categories can be derived from representa¬ 
tions of the quantum group Uq (su(2)) (called the quan¬ 
tum universal enveloping algebra of su(2)). When q 
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is taken to be a primitive root of unity, C/q(su(2)) has 
finitely many irreducible representations with nonzero 
quantum dimensions, which lead to symmetric 6j- 
symbols. An efficient way to construct this data is 
through the Jone-Wentzl projectors in Temperley-Lieb 
algebra (see ref^ for example). Examples include semion, 
Fibonacci, and Ising data, which we discuss now. 


The four quantum double labels are 00, 02, 20, 22. The 
S matrix is 


^1 (j) (j) \ 

(j) —1 (j)^ —(j) 

4> (jP —1 —<j) 

\(jp -<t) -<t) 1 / 


(119) 


1. Doubled semion model 


The twists for the labels 00,02,20 and 22 are 
1, exp(47ri/5), exp(—47ri/5) and 1, respectively. 


Semion data can be obtained at the g-deformation pa¬ 
rameter q = exp(7ri/3). String types are L = {0,1} 
(sometimes denoted by ( 1 , sj), with quantum dimensions 
do = 1 and di = — 1. It has the same fusion rule (5iio = 1 
as that of the group Z 2 representation theory. 

The nonzero symmetric 6j-symbols are 


^000 

^000 


1,G' 


000 

111 


= i,G 


oil _ 
oil — 


(115) 


3. Doubled Ising model 

Ising data can be obtained at the g-deformation pa¬ 
rameter q = exp(37ri/4). 

The string types are L = (0,1, 2}, sometimes denoted 
by {l,cr, "0}. The quantum dimensions are do = l,di = 
\/2 and d 2 = 1. 

The fusion rules are 


The other nonzero 6j-symbols are obtained through the 
tetrahedral symmetry. 

The nontrivial R matrix is R^i = i. 

There are four quantum double labels: 00,01,10,11, 
called boson, semion, anti-semion, and doubled semion. 
The S matrix is 


/ 1 -1 -1 l\ 
- 1-1 1 1 
-11-11 
Will/ 


(116) 


<^000 — l,<5oil — 1,<5 o22 — 1,^112 — 1, (120) 


and the nonzero 6}-symbols G are given by 


>^000 _ 1 ^000 

^000 ~ ^111 


^ ^000 _ 1 _ 


1 


^011 _ ^ ^011 

'-'122 — ^/2 ’ '-'211 


^022 

v/2’^022 


1,G 


112 _ 
112 — 


1 

71' 


( 121 ) 


The other nonzero symmetrized 6j-symbols are ob¬ 
tained through the tetrahedral symmetry. The nontrivial 
R matrices are 


The twists are 1, i, —i, 1. 


2. Double Fibonaeci model 


Fibonacci data can be obtained at the g-deformation 
parameter g = — exp(7ri/5). 

The string types are L = (0, 2}, sometimes denoted by 
{1,t}. Let (j) = 1-^^ be the golden ratio. The quantum 
dimensions of 0, 2 are do = 1 and d 2 = </>. 

The fusion rules are 


^000 — ^022 — ^222 — 1 , ^002 — 0 


(117) 


and the nonzero 6j-symbols G are given by 


^000 _ 1 ^022 _ ^022 

^000 “ ^022 “ ^222 


!/</>, 


^000 _ 

^222 — 


l/v^,Gi2 = 


-I/ 4 P. 


(118) 


The other nonzero symmetrized 6j-symbols are ob¬ 
tained through the tetrahedral symmetry. The non¬ 
trivial R matrices are i ?22 = exp(—47rz/5) and i ?|2 = 
exp(37ri/5). 


i ?22 = “Ij-R^i = “b-^ii = exp(—7ri/8), = exp(37ri/8). 

1122) 

There are nine quantum double labels: 00, 01, 02, 10, 
11, 12, 20, 21, and 22. The S matrix is 


1 

72 

1 

72 

2 

72 

1 

72 

1 \ 

72 

0 

-72 

2 

0 

-2 

72 

0 -72 

1 

-72 

1 

72 

-2 

72 

1 

-72 1 

72 

2 

72 

0 

0 

0 

-72 

-2 -72 

2 

0 

-2 

0 

0 

0 

-2 

0 

2 

72 

-2 

72 

0 

0 

0 

-72 

2 -72 

1 

72 

1 

-72 

-2 

-72 

1 

72 1 

72 

0 

-72 

-2 

0 

2 

72 

0 -72 

1 

-72 

1 

-72 

2 

-72 

1 

-72 1 / 

(123) 


The twist is determined by R matrix. 


X. ELECTRIC-MAGNETIC DUALITY IN 
TOPOLOGICAL THEORY WITH FINITE 
GAUGE GROUPS 


Here we give a consequence of our results in the context 
of electric-magnetic duality in topological gauge theory 
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with a finite (gauge) group H, which has been proposed 
and studied i n^^’^° . 

For an arbitrary finite group H, there is a procedure^^ 
to derive its unitary symmetric 6j-symbols equipped with 
an R matrix, which can be used as the LW input data to 
construct the Repn model, see Sec. IIX Al On the other 
hand, we can define a Vecu model on the same trivalent 
graph, denoted by F. This Vecu can be identified with 
Kitaev’s quantum double model with H defined on a tri¬ 
angulation - the dual graph of F. As discussed in the 
previous section, in the Repn model, fluxons at plaque- 
ttes of F are labeled by conjugacy classes A oi H, and 
charges at vertices (of F) by irreducible representations 
of H. In the Veen model, charges at triangular plaque- 
ttes of the triangulation are labeled by conjugacy classes 
A, while fluxons at vertices (of the triangulation) by irre¬ 
ducible representations of H. In terms of local operators. 
Bp (and Qv) in the Repn model can be identified with 
Qy (and Bp, rep.) in the Vecn model. 

This gives rise to an electric-magnetic transforma¬ 
tion (EMT) between these two models^^. The electric- 
magnetic duality asserts that the two models connected 
by the EMT are actually equivalent to each other—. 

Since the existence of a transformation in general does 
not imply the existence of a corresponding symmetry or 
invariance, the validity of EMD is much stronger than 
the existence of EMT. Well-known examples in quan¬ 
tum field theory include spontaneous symmetry break¬ 
ing and non-abelian gauge anomaly in quantized chiral 
gauge theory^. Namely, one needs to check that the 
EMD in topological gauge theory is not violated by sym¬ 
metry breaking or global excitations. Even if sometimes 
the arguments for the EMD is intuitively simple, the con¬ 
crete checks for exact duality may be highly nontrivial. 
Here we provide two concrete checks for the EMD be¬ 
tween the RepH model and the Vecn model. 

Our first check is to verify that the Hilbert space of the 
two models connected by the EMT has the same dimen¬ 
sion. Certainly this is a necessary condition for the two 
models to be equivalent to each other. We slightly extend 
the RepH model by enriching its Hilbert space again at 
each tail. To each vertex v we associate a tail labeled 
by Qy and a matrix index of representation qy , which 
take values 1,2,... ,dim,j^. (Recall dim(g^) = aq^dq„.) 
Define the Hamiltonian by 

H = -Y,Qy-Y,Bp, (124) 

V p 

where Bp is defined in Eq. (l68l) . The operators Qy and 
B will not affect niy. We still call this slightly extended 
model a Repn model. The Hilbert space is illustrated in 
Eig. [nKa). 

The VecH model is identified with Kitaev’s quantum 
double model with H. The Hilbert space is spanned by 
group elements in H at all edges in the triangulation. See 
Fig.IIZKb). 

The Hilbert space of two models has the same dimen¬ 
sion. To see this, we look at local Hilbert space TLy at 




FIG. 17: (a). The Hilbert space for Repn model on trivalent 
graph r, with each tail labeled by , rriv. (b) The Hilbert 
space for V cch model on the triangulation (solid line) dual 
to r, with local Hilbert space C[H] on each edge. 


each vertex. It has basis {\i, j,k,l,q,m)}, labeling the 
following diagram 


I gm 



where i,j,k are labels on three incoming edges, and 
l,q,m the enriched charge degrees of freedom. Note that 
there is exactly one tail to each vertex. The dimension 
of Ry is 

dim(H„) = SjkApq.dq = didjdk- (125) 

Therefore, effectively the local Hilbert space at each edge 
e labeled by je has dimension d^^. By a theorem for the 
order (number of elements, or dimension of the group 
algebra) of a finite group, the local Hilbert space at each 
edge has dimension This gives the same 

dimension of the Hilbert space of the Vecn model on the 
triangulation. In fact, there is a duality transformation 
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between the two models 
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with the summation {j} over labels on all edges on 
graph r (or on the triangulation), and Vj for represen¬ 
tation space of j. The map on the 4th line is an iso¬ 
morphism composing two Wigner’s 3j-symbols or the 
Clebsch-Gordan coefficients that decomposes the tensor 
product of two representations. The map on the last 
line is a generalized Fourier transformation (i.e. the 
Peter-Weyl theorem) between two bases of C[H], from 

to \^)h^H- 


Paj3 




Our second check is about the spectrum of the models. 
The equality of the ground state degeneracy in the Repn 
and VecH models has been verified in, e.g.^. With our 
quantum double classification of quasiparticle excitations 
in the Repn model, we are able to check the equivalence 
of the two models at the level of quasiparticle excita¬ 
tion species. Previously it is known that in the Vecu 
model (Kitaev’s quantum double model) it is the quan¬ 
tum double of the finite group H that classifies the ele¬ 
mentary excitations^. In this paper we have shown that 
the RepH model (or the LW model with input data from 
finite group H) accommodates dyon excitations classified 
by the quantum double of H. So the two models have 
the same excitation spectrum in their quantum numbers 
and energy levels. This is certainly a highly nontrivial 
check for the electric-magnetic duality between the two 
models. 


XI. RELATION TO TOPOLOGICAL QUANTUM 
FIELD THEORY 


A. GSD and Turaev-Viro TQFT 

We first consider zero temperature case. We denote the 
input data {d, S, G} by a unitary fusion category C that 
derives them. The zero temperature partition function of 
Levin-Wen models on surface S equals the ground state 
degeneracy (GSD). The GSD is is related to Turaev-Viro 
invariant - a topological invariant for 3-manifold defined 
below - by 


GSDc(S) = X S'!). (127) 

We first define and then sketch the proof. 

Given C and a compact oriented 3-manifold M, we con¬ 
struct the number {M) as follows. Any 3-manifold M 
has a triangulation, i.e., can be discretized into tetrahe¬ 
dral. We choose an arbitrary one, and the desired number 
will be triangulation independent. 

1. Assign labels to all edges. 

2. Assign 6j-symbols to all tetrahedral as follows. 
Due to tetrahedral symmetry in Eq. ®, such 6j- 
symbols do not change under rotation of tetrahe¬ 
dral. 


^ (128) 


3. Assign quantum dimensions dj to all edges labeled 
by j- 

4. Assign 1/D to each vertex (in the triangulation). 

5. Multiply all the quantities in step 2-4, and take 
the product over tetrahedral, edges, and vertices of 
these numbers. 

6. Sum over all labels. 

We get 

rPiM) = ^ n n d n (6j-symbols). 

labels vertices edges tetrehedra 

(129) 

This number does not depend on choices of triangula¬ 
tion. 

Now we sketch the proof in Eq. (I127D . The ground 
states are Hp Bp = 1 eigenvectors. Hence GSD = 
trdlpl^p)- To relate the trace to tJ^, we first write 
down the dual triangulation of the trivalent graph as fol¬ 
lows. 



Levin-Wen model is viewed as a Hamiltonian ap¬ 
proach to Turaev-Viro topological quantum field theory 
(TQFT). The topological observables in the former are 
related to topological invariant of 3-manifolds. 



( 130 ) 
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Eq. m becomes 





( 131 ) 


It can presented by three tetrahedral as follows. The top 
three triangles are those in the bra of above equation and 
the three bottom ones are in the ket. 



This identihcation with leads to tr(J|pi3p = 


B. Excitations and the extended Turaev-Viro 
invariant 


In this subsection we will explain how the excitation 
in the above model are related to an extension of the 
Turaev-Viro invariant to manifolds containing links de¬ 
fined by Turaev and Virelizier— . 

Here we continue using the conventions of the last sub¬ 
section. Let Z{C) be the quantum double category asso¬ 
ciated to C. As mentioned above a minimal solution to 
Equation (1351) is identihed with a quantum double ele¬ 
ment J. Let L be a link in E x whose components are 
labeled with such quantum double elements. 

Iii2^, the Turaev-Viro invariant is extended to Z{C)- 
colored links and in particular, defines an invariant of 
the pair (T, x S^, L). In this work the invariant is defined 
using skeleton which are 2-polyhedron with certain prop¬ 
erties. Taking the dual of a triangulation gives a skeleton. 
A link in a skeleton is a collection of loops immersed in 
the 2-dimensional simplices of the skeleton with certain 
transversality conditions. A quantum double label J and 
its associated half braiding z'^ can be used to define a new 
symbol similar to a 6j-symbol, see^. The extended T-V 
invariant is defined in a similar way to the T-V invariant 
outlined above: First, a skeleton of (E x S^,L) can be 
decomposed into building blocks (which are analogous to 
tetrahedron). Each face is assigned a label (here the a 
face of the skeleton is dual to an edge in the triangu¬ 
lation). Each building block correspond to a 6j-symbol 
or a new symbol coming from the quantum double la¬ 
bels of the link. As in Eq. (11291) . the extended invariant 
is obtained by taking a weighted sum over all possible 
labelings of the product of these symbols. 

The model given in this paper is a Hamiltonian realiza¬ 
tion of the extended T-V invariant. As outlined below, 
the link L in T, x can be associated with an operator 


/l which is a composition of certain operators given in 
Section |Vl (Note /l is not unique.) 

For simplicity, let us first describe the situation when 
E = 52 . In Fig. HlKa), the bottom plate presents an 
initial state with two dots presenting dyons. The top 
plate presents a final state. The operator is defined 
by composing certain creation and annihilation opera¬ 
tors and charge contractions, with the composition order 
coinciding with the time direction in the figure and de¬ 
termined by the topology of L. Since the string operator 
is path independent, only depends on the topology of 
the tangle, i.e. the portion of the link between the two 
plates. The operator /l is parameterized by charges at 
the ends of strings on top and bottom plates. 

We require that the closer along of the tangle un¬ 
derlying /l is the link L, see Fig. fl^ cL When the link 
is contained in a 3-ball the closer is trivial and the tangle 
can be chosen to be the link. For example, when L is the 
Hopf link Hj k in a 3-ball in E x whose components 
are labeled with J and K then is the composition of 
the operators described in Figure [8l 

In general, E and L may not be in a 3-ball. 

The description of /l in such a situation involves the 
topology of E and/or a nontrivial closer of the tangle 
along S^. However, since the operators are all local, 
can be described by a similar process as above. 

The extended Turaev-Viro invariant associated to 
C is equal to the Reshitkhin-Turaev invariant ^ 

S^,L) associated to the quantum double category Z{C), 
see^. The R-T invariant is defined by representing the 
3-manifold by surgery on some framed link K in then 
applying certain quantum invariants coming from Z{C) 
to the link K U L. Thus in the definition of the R-T 
invariant one does not have to work with trivalent graphs 
but the evaluation of certain link invariants. 

An argument similar to the one in the last section 
shows that 

tr(/L) = rc^'^(E X S\L) = x S\L). (132) 

The above trace is taken in the Hilbert space and is the 
sum over all charges on the open strings. This trace can 
be viewed as a charge contraction in the time direction, 
and connects all open strings in Fig. fTSI a') to a closed 
link in Fig. [THT cL In the case of ground states i.e. when 
L is trivial Eq. (jl32l) is just Eq. (I127D . 

In particular, when the link L is trivial, /l = r ip 
is the ground state projection operator, and Eq. (11321) 
recovers Eq. (11271) as a special example. 
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(b) (c) 


FIG. 18: (a) A tangle representing /l, which is the creation 
of a dyon-pair from a ground state, (b) Illustration of gluing 
dots along S^. (c) Closure of underlying topological object is 
L. 


XII. CONCLUSIONS AND DISCUSSIONS 

In this paper we have studied how to describe the full 
spectrum of dyon excitations in the generalized Levin- 
Wen models. Previously it was known that in the LW 
models, fusion of two pure fluxons generally may lead to 
the appearance of charge d.o.f. To incorporate the latter 
explicitly, we enlarge the Hilbert space by introducing a 
tail (labeled by a string type) at one of the edges of each 
vertex, and modified the LW Hamiltonian accordingly. 
Though we have to deal with new configurations with an 
extra tail at each vertex, in this approach we have been 
able to achieve the following: 

(1) In our extended Hilbert space with enriched d.o.f. 
for charge at vertices, we are able to study the proper¬ 
ties of charge and fluxon-type of dyon excitations, and 
in particular their interplay through the twist operation. 
We have shown that one needs three quantum numbers 
“ charge, fluxon-type and twist - to describe the dyon 
species for elementary excitations, or the total dyonic 
quantum numbers of excitations localized in a region 
“surrounded by vacuum”. We emphasize the necessity 
of introducing the twist, as the third quantum number 
beyond the charge and fluxon-type, for a complete de¬ 
scription of a dyon species. 

(2) The above conclusions are obtained by studying 
the operator algebra formed by local operators and its 
irreducible representations (simple modules). We have 


shown that all local plaquette operators preserving topo¬ 
logical symmetry, i.e. invariant under Pachner moves 
form the so-called Tube algebra. The latter is a general¬ 
ization of usual Bp operators in the LW models; in fact, 
the operators Bp form a subalgebra of the Tube algebra. 

(3) String operators can be realized as linear maps on 
the extended Hilbert space. Irreducible representations 
(simple modules) of the Tube algebra are shown to be 
in one-to-one correspondence with the half-braiding ten¬ 
sors that are used to define string operatorsi2i2^. In this 
way, we establish that the Tube algebra and string op¬ 
erators are dual to each other by a (generalized) Fourier 
transformation. On the other hand, half-braiding tensors 
are ingredients to define the quantum double (the center) 
category of the input unitary fusion category. So we con¬ 
clude that the quantum numbers of dyon excitations are 
organized by irreducible representations of the Tube al¬ 
gebra, or equivalently by the quantum double category, 
as the center of the input unitary fusion category of the 
LW model. Twist is a property necessarily associated 
with the quantum double category. 

(4) Realizing string operators, as linear maps over the 
extended Hilbert space enables us to obtain not only the 
S, T matrices, but also the braid group representations 
for dyons. This knowledge is important for describing 
emergent braid statistics of dyon excitations^! and, there¬ 
fore, will play a crucial role in designing quantum com¬ 
putation codes that exploit manipulation of excitations 
in the topological phases. 

(5) We can systematically construct explicit 
states/wavefunctions with given quantum numbers. 
This enables one to study more physically interesting 
quantities, such as entropy, entanglement entropy etc, 
and to design quantum computation algorithm based on 
manipulation of the non-abelian anyonic quasiparticles. 

(6) A consequence of our results is that the Kitaev 
quantum double model (the Toric Code model) associ¬ 
ated with a finite group on a triangulation and the (gen¬ 
eralized) LW model with input data from the same fi¬ 
nite g roup has the exactly same dyon excitation spec¬ 
trum, characterized by the same quantum double cat¬ 
egory. This provides a strong check/test/eveidence for 
the electric-magnetic duality between the two models, 
not only for ground states but also at the level of the full 
excitation spectrum. 

Finally we want to emphasize the following point. One 
may add more terms into the LW Hamiltonian, which 
may not commute with the existing two terms. When 
the coupling strengths of these additional terms are suf¬ 
ficiently small, we expect that the model remains in 
the same topological phase, with the energy levels of 
the many-body states getting shifted, provided there is 
no level crossing between the ground states and excited 
states. With such more general Hamiltonians, we be¬ 
lieve that the quantum double category or the pertinent 
Tube algebra of local operators we have obtained for the 
LW model still provides a “complete basis” for many- 
body excitation states in the enlarged Hilbert space and. 
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therefore, could still be useful. For example, we may use 
this “basis” to formulate/compute perturbation theory 
corrections for elementary excitations. 
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Appendix A: Some properties of quantum double 

The quantum double category is characterized by the 
half-braiding tensors z. We list some properties and sym¬ 


metry conditions on z. 

Orthonormal relation: 

^fjqt ^Ijpt ~ ^pq^jpt* , 
i 

^qjlt ^pjlt ~ ^pq^pjt* , (Al) 

I 


satisfies the symmetry conditions 


Y - \ ' d G-TYr" J 

'pjqt ^gj*pr 

r 

(A2) 

,J — N ' d z'^ 

'qj-pr ^ -‘pjijt J 

(A3) 


t 


where the second condition is a consequence of the first 
one together wit the orthogonality relation ([6|). 
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